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SUMARY 


A simplified method has been devised for the solution of two- 
dimensional compressible flows through well-defined passages. This 
method makes use of plastic cams which automatically set the length- 
to-width ratio of rectangles formed by streamlines and equipotential 
lines represented by spring-steel wires. Pressure distributions around 
four cascades of turbine blades and along the surface of a choked nozzle 
determined by this method are shown to compare well with experimental 
results. 


INTRODUCTION 


Two-dimensional flow problems may be solved by the potential -flow 
plotting method. Reference 1 describes a method of determining the 
velocity distribution for a cascade of turbine blades for incompressible, 
inviscid flow by the use of a wire-mesh flow plotting device. Since the 
local velocities on gas-turbine blades commonly reach sonic velocity, 
the effects of compressibility should be taken into account when solving 
for the blade-surface velocity distribution. The problem of plotting 
compressible flow is more difficult than that of plotting incompressible 
flow since the density is not constant but varies with velocity. For 
that reason, the rectangles formed by streamlines and equipotential lines 
of a flow plot are not curvilinear squares as in the incompressible case 
but vary in length -to -width ratio as the velocity varies in the field. 

In previous attempts to plot compressible flows, the length-to-width 
ratio of the rectangles has been adjusted individually. Brenner and 
Kilgore, in reference 2, attempted a solution by manually measuring the 
width and setting the length of each rectangle. Sells, of General Electric 
(reference 3)> devised four -prong calipers which set the correct length 
when adjusted to width. In these methods the streamlines were adjusted 
individually in a time-consuming iterative process. 

This paper describes a cam device which continuously maintains the 
correct length-to-width ratio of each rectangle formed by the streamlines 
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and equipotential lines which in this case are represented by spring-steel 
wires. The compressible, potential flow through well-defined passages 
such as high-solidity cascades, nozzles, and ducting can he obtained 
directly with this device. Comparisons of surface pressure distributions 
obtained with this device with high-speed cascade data are shown for 
four turbine -blade cascades. A two-dimensional nozzle shape redesigned 
with this device in order to avoid supersonic velocities along the sur- 
face of the nozzle was also studied. 


SYMBOLS 


a 

C 

M 

P 

q 

V 

X 

y 

AX 

AY 

P 

7 

0 

<t> 

P 


speed of sound 
constant 

Mach number (V/a) 
static pressure 

dynamic pressure G pv 2 ) 

velocity 

distance parallel to flow direction 
distance normal to flow direction 
distance between equipotential lines 
distance between streamlines 

inlet flow-angle, angle between upstream velocity and 
normaVto stagger line 

ratio of specific heats 

turning angle, angle between entering velocity and 
exiting velocity 

potential function 

density 
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Subscripts: 


1 


upstream 


2 


downstream 


o 


stagnation 


cr 


critical value; that is, value when local Mach number 
is 1.0 


PRINCIPLES OF COMPRESSIBLE -FLOW PLOTTING 


Design and construction of device .- In order to satisfy the equa- 
tions of compressible potential flow, the flow plot must satisfy bound- 
ary conditions, have the length-to -width ratio of all rectangles equal 
to the ratio of local density to stagnation density, and maintain orthog- 
onality of all intersecting streamlines and equipotential lines. 

The compressible -flow plotting device consisted of a grid of wires 
representing streamlines and equipotential lines, pins at the inter- 
section of the wires, and plastic cams which restrain the pins so that 
the rectangles formed by the wires have the proper length-to -width 
ratio. Figure 1 shows the arrangement of the cams, pins, and wires in 
a small element of the compressible-flow plotting device. 

The wires were spring -tempered stainless steel. A diameter of 
0.04l inch was selected for the wires as a compromise between the flex- 
ibility needed to meet boundary conditions and the stiffness needed to 
keep the streamlines and equipotential lines smoothly faired. The pins 

were of — -inch-diameter brass rod, inches in length. Two holes were 

8 ^ 

drilled in each pin for the wires 1/8 inch and 1/4 inch from the base 
at right angles to each other. 

The cams (fig. 2) were manufactured by injection molding of a 
phenolic -resin plastic. The thickness of the cams was 0.060 inch and 

the width of the cam slot for the pins was 0.125 t oiooi ^ nc ^ 1, The 

coordinates were accurate to ±0.003 inch. The range of Mach number for 
the cams was from 0.3 to 1.0. The maximum width of streamtubes for 
which the cams were designed was I .760 inches and the minimum width was 
0.864 inch. Small pointed pins were strung along the boundary stream- 
lines and were pressed into a soft-wood board to fix the boundary con- 
ditions. 
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The shape of the slots of these cams was so designed that the 
width AY is the proper function of length AX. Since d$ = V dx 
and the change in potential between adjacent lines of equal potential 
is constant, then 


V 


_C_ 

AX 


( 1 ) 


where C is an arbitrary constant chosen to give a convenient cam size. 
As already noted. 



( 2 ) 


For isentropic flow 


£_ 



(3) 


Substituting equations (l) and ( 3 ) into equation (2) and solving 
for AY gives 

AY = ^ (A) 



The relation between AX and Mach number may be determined from the 
energy equation 


|V2 + 




or, after the terms are rearranged. 



(5) 
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By substituting equation (l) into equation (5), the Mach number in terms 
of AX is 


M 


2 



( 6 ) 


Equations (4) and (6) are plotted in figure 3 for a value of C = — 

2 

which gives a cam of convenient size if AX is measured in inches. As 
can be seen from this figure, the values of AX are difficult to deter- 
mine accurately for values of AX near M = 1.0. The cam slots were 
designed from equation (4) for 7 = 1.40. The coordinates measured from 
the cam center for each slot are AX/2 and AX/2 (fig. 2)> hence, the 
distance between identical points on the four slots is AX and AX. Each 
value of AX and AX on the cam corresponds to one Mach number which is 
given by equation (6). The cam slot was not designed to include supersonic 
velocities because the discontinuities present in supersonic flow cannot 
be treated. 

Although it has not been shown mathematically, it was found experi- 
mentally that orthogonality was maintained by the cams alone with two 
exceptions, which are shown in figures 4(a) and 4(b). Although no wires 
were used with the pins and cams to obtain these patterns, lines connecting 
the intersections have been drawn in the figures to show the pattern of 
streamlines and equipotential lines. The first figure shows distortion 
which is equivalent to a displacement of the streamlines along a 45° line. 
In the second pattern (fig. 4(b)) the cams turned alternately and gave 
zigzag streamlines. Adding the wires prevents these extraneous patterns 
from occurring by smoothly fairing the streamlines and equipotential 
lines. The wires are also used to keep the pins vertical. 

Since the cams are finite rather than infinitesimal in size, a 
difference exists between the distance measured along the curved stream- 
line and the straight-line distance between pins. This difference depends 
on the radius of curvature and is insignificant for these cams for radii 
of curvature greater than 3 inches. Hence, any plot made must be of 
sufficient size so that a small radius of curvature is not present on any 
important boundary line and enough streamtubes are included to define the 
flow throughout the field. 

Procedure . - The compressible-flow plotting device is constructed on 
a scale drawing of two passages of a cascade. One passage is sufficient 
if the boundary streamlines are made identical, but it is better to plot 
the flow in two or three passages to a point far enough past the nose to 
insure that the flow around the nose of the blades is correctly defined. 
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The scale of the drawing is uniquely determined by the scale of the cams, 
the number of streamtubes per passage, and the entering Mach number. In 
general, the scale of the drawing must be changed in order to vary the 
entering Mach number. The distance between two adjacent stagnation stream- 
lines must be equal to the value of AY for the entering Mach number times 
the number of streamtubes between stagnation streamlines. Thus, by 
determining the distance between stagnation streamlines upstream, the 
scale of the drawing is found. 

The complete grid of cams cannot be easily adjusted because of the 
great amount of friction on the cam surfaces. Hence, for each compressible- 
flow plot to be made, the grid of wires and pins should be shaped to the 
boundary conditions before adding the cams. The cams are added in rows 
starting from the uniform upstream flow. This procedure eliminates 
interference between cams since each cam is on a different level than 
the adjacent ones. Because of this procedure, the cams are slightly 
tipped both in the stream direction and normal to it (fig. 1). Measure- 
ments are made of the values of AX along the boundary lines and the 
Mach number is found for each measurement by equation (6). From a measured 
value of AX an average value of the velocity for a finite distance 
along the boundary can be found. This average velocity is considered to 
occur midway between the pins. 

A check on the final result is that the normal-force coefficient 
computed from the integrated pressure distribution around the blade must 
correspond to the turning of the air through the cascade. Actually, 
the amount of turning through a given cascade is, to some extent, arbitrary; 
that is, in general, the inlet air angle, inlet Mach number, and cascade 
configuration are not sufficient to determine the exit air angle. Spec- 
ifying the position of the trailing-edge stagnation point would make the 
problem unique. Assuming a local flow direction at the trailing-edge 
region is equivalent to assuming a position of the trailing-edge stag- 
nation point. When the compressible -flow plotting device is used, the 
existence of boundary layer is recognized and the exiting flow is assumed 
to have no velocity peaks around the trailing edge; hence, the local 
streamlines leave with approximately the mean-line direction and with a 
wake thickness equal to or greater than the trailing-edge thickness. In 
the solutions that have been obtained with this device, this direction 
has been assumed as the exit direction for the first approximation, with 
some subsequent small adjustment of this direction, if necessary, in order 
to satisfy the criterion just mentioned. The adjustment is readily 
accomplished since it affects the turning angle directly but has almost 
no effect on the normal-force coefficient as computed from the pressures. 

For high-solidity cascades with relatively low entering velocities, the 
position of the nose stagnation point has only a small effect on the 
pressure distribution. For this reason, no attempt was made to determine 
the local velocities at the nose accurately. 
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COMPARISON WITH EXPERIMENTAL RESULTS 


The device has been used to plot the flow through four turbine-blade 
cascades , and the pressure distributions measured from these plots are 
compared with those from high-speed cascade tests. Figure 5 shows the 
pressures obtained along the surface of an airfoil in an 88° turning 
cascade for an entering Mach number of approximately 0.34 for both the 
test and the flow plot. It can be noted from this figure that an irreg- 
ularity of points occurs in the region of high velocity. This irregu- 
larity can be attributed to the large changes in AX (or velocity) that 
occur for small changes in AX at speeds approaching a Mach number of 1. 

A photograph of the plotting device for this cascade is shown as figure 6. 
Several cams were left out along the pressure surface of the blade because 
the Mach number in that region was below 0.3, which is the minimum value 
for which the cams were designed. In small low-speed areas such as this, 
the length -to -width ratio of the rectangles is adjusted manually rather 
than by using special large-size cams. The rectangles will become almost 
square as in an incompressible flow with a density ratio of 1.0 when the 
flow velocity is very low. 

Figure 7 compares the pressure distributions obtained from a plot 
and from the cascade test for the same blade section as in figure 5 but 
at a higher inlet air angle. The difficulty of maintaining the boundary 
conditions around the highly curved nose and the high velocity in this 
region caused the difference in the pressure distributions. 

Pressure distributions for a cascade similar to that shown in figure 5 
are compared in figure 8. The same small amount of irregularity occurred 
in the pressure distribution obtained by the flow plotting device shown in 
this figure. There is no irregularity of points in the pressure distribu- 
tion for the fourth cascade (fig. 9) which was a low -turning turbine 
cascade. Evidence in other plots also has shown that this irregularity 
of points occurred only in regions of high velocity and did not occur 
in a steep velocity gradient. Generally good agreement is shown between 
the pressure distributions of the cascade tests and the flow plots for 
the four cascades tested in spite of the approximations made at the 
nose. 


The average length of time required to plot the flow through a 
turbine cascade is 8 hours, an appreciable reduction in time from that 
required with the incompressible -flow plotting device. The time reduction 
is due to the setting of the streamlines and equipotential lines by the 
cams. The accuracy of the plots made with the compressible -flow plotting 
device is generally good, results comparing equally as well as those of 
the incompressible-flow plotting device with cascade data. Reliability, 
of the results should be much better for the cam device since the plot 
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required no alinement, positions of the streamlines being rigidly fixed 
by the cams. 

All the plots were made for entering Mach numbers which were slightly 
below the critical value for the blade section so that the maximum Mach 
number along the blade surface would be less than 1.0 since the mechanical 
inaccuracies of the cams produce larger errors near sonic velocities. 

An effort has been made to determine the entering Mach number for which 
sonic velocity will first occur at some point along the blade surface. 

This Mach number is of value in comparing various blade sections. 

A simple method which has been used to obtain the critical entering 
Mach number is to extrapolate the highest Mach number obtained on the 
blade surface in the flow plot to a Mach number of 1.0 by means of the 
Prandtl-Glauert extrapolation method (reference 4) . To use this method, 
a mean passage Mach number is first calculated for the section of the 
passage where the maximum surface velocity occurred in the flow plot. 

This mean passage Mach number is best obtained from the upstream Mach num- 
ber and the ratio of the area of flow upstream to the area of the passage 
at the point of maximum local Mach number. The mean passage Mach number 
is considered to be the equivalent of M^ in the Prandtl-Glauert 
extrapolation which is equation (8-4l) of reference k 


C 


P 



p - p 

where C-n = and Cr> is the low-speed or incompressible value 

P loo 

of the pressure coefficient. The extrapolation is made for a constant 

value of C-n from the condition of the flow plot to (cJS which is 

Po \ P/ cnt 

given by equation (3-21) of reference 4 


( c p) 


crit 


7M 0 


+ 1 


tt"-; 


7 

2 \~- T 


- 1 


A single plot of these two equations facilitates their rapid simultaneous 
solution for the critical value of Moo or critical mean passage Mach 
number. The value of the mean passage Mach number at this critical 
point is then converted to a critical entering Mach number by using 
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the original area ratio. Critical entering Mach numbers computed by 
this method for the four flow plots and critical Mach numbers obtained 
in cascade tests are as follows: 



Mcr 

First 

configu- 

ration 

Second 

configu- 

ration 

Third 

configu- 

ration 

Fourth 

configu- 

ration 

Cascade test 

0-354 

0.385 

0.467 

0.473 

Compressible -flow plot 

.356 

.395 

.524 

.46o 

Maximum blade -surface 
Mach number 

.838 

.848 

.867 

.850 


Very good agreement of the critical values is obtained except in the 
case where peak velocity obtained with the plotting device was in error 
(third configuration). 

The compressible -flow plotting device was also used to determine 
the velocity distribution along the surface of a choked subsonic 
two-dimensional nozzle. The solution could not be made for the complete 
nozzle shape since the boundary streamline could not be made to conform 
to the nozzle surface where it was highly curved even though the local 
Mach number was 1.0. This deviation indicates that the velocity over the 
nozzle surface would have been supersonic. Hence, the boundary stream- 
line along the surface was allowed to deviate slightly from the nozzle 
in the supersonic region to increase the radius of curvature so that 
supersonic velocities were avoided (fig. 10). Subsequent test of the 
nozzle shape confirmed the existence of a supersonic region along the 
surface of the nozzle where the compressible -flow plotting device indi- 
cated supersonic velocities. Surface Mach numbers along the unmodified 
section of thje nozzle show reasonable agreement considering the fact 
that the velocities were very close to sonic and the compressible-flow 
plot was made without consideration of the boundary layer. 


CONCLUDING REMARKS 


The plotting device described herein simplifies the graphical solution 
of two-dimensional compressible-flow problems that would otherwise require 
long and tedious mathematical solutions. The cams used on the device 
greatly speed the solution by positioning the streamlines and equipotential 
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lines of the flow plot and compensating automatically and continuously 
for the change in density throughout the flow field. 

Reasonable results were obtained for four turbine-blade cascades 
and a two-dimensional nozzle. Other configurations of the channel -flow 
type such as ducting, turning vanes, and passage shapes which are long 
compared to their width could be investigated by use of the compressible - 
flow plotting device. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., January 15, 1952 


REFERENCES 


1. Westphal, Willard R., and Dunavant, James C.: Application of the 

Wire-Mesh Plotting Device to Incompressible Cascade Flows. 

NACA TN 2095, 1950. 

2. Brenner, Wm. C., and Kilgore, Lee A.: Air Flow Determination by 

Graphical Methods. Preprint, Inst. .Aero. Sci., Inc., 1944. 

3. Poritsky, H., Sells, B., and Danforth, C. E.: Graphical, Mechanical 

and Electrical Aids for Compressible Fluid Flow. Data Folder 46150, 
Aircraft Gas Turbine Eng. Div., Gen. Elec. Co., Aug. 13, 19^* 

4. Liepmann, Hans Wolfgang, and Puckett, Allen E.: Introduction to 

Aerodynamics of a Compressible Fluid. John Wiley & Sons, Inc., 

19^7, PP- 28, 244-245. 



Figure 1.- A small element of the compressible -flow plotting device. 
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2.- Plastic cam and ordinates. 


AY, in 


NACA TN 2681 


13 


I 

I 



AX, in. 


c_ _ 1 

a o ” 2* 


Figure 3*- Mach number for distances on cam. 
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(a) Displacement along 45° line. 


Figure 4. - Exceptions to orthogonality. 
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(b) Zigzag of streamlines. 
Figure 4.- Concluded. 
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Figure 5«- Pressure distribution for turbine blade. 
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Figure 6.- Compressible-flow plot of turbine blade of figure 5* 
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Figure 7.- Pressure distribution for turbine blade. 
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Figure 8.- Pressure distribution for turbine blade. 
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Figure 9.- Pressure distribution for turbine blade. 
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Figure 10.- Nozzle shape and Mach number distribution along nozzle surface 

in choked condition. 
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